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Abstract 



The Bose gas in an external potential is studied by means of the local density 
approximation. An analytical result is derived for the dependence of the critical 
temperature of Bose-Einstein condensation on the mutual interaction in a generic 
Tj- ' power-law potential. 
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In an ideal Bose-Einstein gas, the zero momentum state can become macroscopically occupied, 
^\ [ and the system then undergoes a phase transition — the Bose-Einstein condensation (BEC) Q. 

The A-transition of liquid helium-II is often interpreted essentially as a result of BEC, the strongest 
degeneracy effect of a Boson system[p. For many years, it was considered hopeless to experimentally 
observe BEC in an atomic gas with weak interactions. With the development of techniques to 
' trap and cool atoms, BEC was recently observed directly in dilute atomic vaporsj^, ^. The new 
^ ' experimental achievements have stimulated great interest in the theoretical study of inhomogeneous 
Q , Bose gases. 

For BEC in the atom-trapping experiments, the critical temperature Tc can be altered by the 
spatially varying potential of trapping. Furthermore, the noninteracting Bose gas can lead to 
unphysical results. The interaction between atoms may have a significant effect on the critical 
I temperature. There have been several investigations analyzing the dependence of the critical tem- 
perature on the trapping potential and weak interaction in the Bose gas[^, 0. Here we shall derive 
an analytical result for the critical temperature of a trapped interacting Bose gas by means of the 
local density approximation. 

When the energy level spacing of the trapping potential is much smaller than kT = jS^^ , the 
local density approximation ||8|, |9|, |l^ is a good approximation. Space may then be divided into 
small cells, and in each cell we may consider the trapping potential energy V{r) to be a constant. 
In this approximation for an ideal Bose gas, the local density p(r) of particles in the gaseous phase 
(i.e., without BEC) is 

p(r) = A-3 53/2(0, (1) 

where 
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e = zexp[-/3F(r)], 



9u{x) = "x^, 



and /i is the chemical potential. Equivalently, the density of states at energy e is| 
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J - V{r)d?r. 
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We may now extend the argument of BEC for rigid walls Q to the case of the spatially varying 
trapping potential V{y). From Eq. (2), the total number of particles can be written as 

N = No+ [ p{r)d^r, (5) 
Jv 

where Nq is the number of particles in the ground state. Let us consider a generic power-law 
potential as in Ref. @ 
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Since §3/2(0 ^ bounded, positive, monotonically increasing function of ^, for the potential (^), 
the integral in Eq. (^) reaches its maximum value at fugacity z = 1 or // = . Prom Eqs. (^) and 
(|5|) the critical temperature Tc, at which the ground state begins to take on macroscopic values, is 
determined by 



drg3/2[exp{-PcV{r))] 

exp(-j/3cl^(r)) dr. 



The definition of the Gamma function 



implies 
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For the specific choice (^) of the potential, expression (^) then becomes 

^ 8i'3/2 r(i/p)r(i/or(i/g) 
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where 

i{p,i,q) = (plq)-' r(i/p)r(i//)r(i/g), 

and ((v) = J"'' is the Riemann zeta function and r] = 1/p + l/l + 1/q + 1/2 . Noticing the 
dependence (g) of the thermal wavelength A on the temperature, we finally obtain 
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From Eqs. (|5|) and (10) we can obtain the fraction of condensation at a temperature T below Tc 



No 
N 



1 



(13) 



where A'o is the number of Bose particles in the BEC state. 
For the harmonic potential 

Mt \ ^ 2 2,1 22 

V(y) = -muj^r^ + -miv^r^, 
the critical temperature and fraction of condensation reduce to 
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where we have made use of 



C(3) = 1.202, /(2,2,2)= 7r 



3/2, 
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The results dH) and (||) were first obtained in Ref. |l^ and have been verified in a recent 
experiment ||l2|. 

Relation (|12|) for the generic power-law potential may also be derived in terms of Eq. (^). In 
Ref. p it is given that 



^ i/p i/i i/g 1 



where 
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exp{9) - 1 



de, 



F{pJ,q) 



(1 - xPji/^+i/^+i/' / dx {I - x'y/"+'/^ / dx(l-x«)i/2 
In fact, we may integrate both Q and F to obtain 

F(p,Z,g) =4V^/(p,/,g)/r(7? + l), 
where in deriving the last equation we have used 

It can then be verified that Eq. (^) coincides with (^) as it should. 
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Along similar lines we may analyze the critical temperature of a trapped interacting Bose gas. 
In Ref. [10| it has been given that when the s-wave scattering length \a\ ^ A, in the local density 
approximation the local density of the gaseous phase as given in Eq. (|^) is replaced by 

p{r) = X-'gs/2{i), (24) 

with 

I = zexp[-/5F(r) - 4aX^p{r)], (25) 

where the term containing the scattering length a describes the mutual interaction of atoms. 
From Eqs. and (|25|) we have 

p(0) = A-3g3/2[^exp(-4aA2p(0)]. (26) 

Thus, at the critical temperature the fugacity is 

Zc = exp[4aA^pc(0)] (27) 

where 

Pe(0) = A-353/2(l). (28) 

For this Zc Eq. (p3) becomes 



P(r)Ai? = ff3/2[^cexp(-/3eF(r) -4aA2p(r))] 

^ g3/2[eM-PcV{r))] + Aa~Xl [p,(0) - p(r)] 9i/2[exp(-/3,y(r))] 
~ 53/2 [eM-PcVir))] + 4aA2p,(0) g,/^ [exp(-4F(r))] 

-4aA-i g3/2[exp(-/3eF(r))] <7i/2[exp(-/5eF(r))], (29) 

where we have used ^ 

53/2 (^) = 51/2 (^), (30) 

and kept only the lowest order in a/ A. The determination of critical temperature Tc for the 
interacting Bose gas involves integration of p{r) in (29). After some calculations similar to those 
in the derivation of (jl^) from we find 

iVA^ = 8 [ Civ + 1) + 4aA-i ( CimCiv) " xiv " 1/2) ) ] Jil^^fl/l^lli^, . (31) 

Pc Ci €2 ^3 

where 
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We finally obtain the critical temperature 



'^8m3/2[C(r? + l)+4aAc-^(C(3/2)C(r/)-x(r/-l/2))]LiL2L3/(p,Z,(z), 
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Comparing Eq. ( ^3|) with Eq. (^) for the critical temperature Tc of the trapped ideal Bose gas, we 
can derive the relation between the critical temperatures of BEC for ideal and nonideal Bose gases 
in the power-law potential 



Tc 



C(^ + i) 



(34) 



which leads to 
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Erom this analytical relation between the critical temperature and the scattering length a, we see 
that the critical temperature Tc decreases for a > 0, and increases for a < 0, compared with Tc of 
the ideal Bose gas. This agrees with Ref. Q. 

We may also obtain an expression similar to the case of an ideal Bose gas for the condensate 
fraction at a temperature T < Tc to be 



No / T y+i ^ ^ a C(3/2)C(??) - X{r] - 1/2) f T 



N \TcJ \c C(?? + i) Vr, 

Eor the harmonic potential (|l^), Eq. (|^) reduces to 

~ ^ C(3/2)C(2) - x(2/3) /SfcrnT^'^' 
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which has been obtained in Ref. 10] • 

The main effect of an external potential is to concentrate Bose particles. A repulsive mutual in- 
teraction between particles will weaken the effect of the external potential in concentrating particles 
around the center of the potential, hence lower the critical temperature. 
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